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We study the mechanical buckhng of a two dimensional membrane coated with a thin layer of super- 
fluid. It is seen that a singularity (vortex or anti- vortex defect) in the phase of the quantum order 
parameter, distorts the membrane metric into a negative conical singularity surface, irrespective of 
the defect sign. The defect-curvature coupling and the observed instability is in striking contrast 
with classical elasticity where, the in-plane strain induced by positive (negative) disclinations is 
screened by a corresponding positive (negative) conical singularity surface. Defining a dimension- 
less ratio between superfluid stiffness and membrane bending modulus, we derive conditions under 
which the quantum buckling instability occurs. An ansatz for the resulting shape of the buckled 
membrane is analytically and numerically confirmed. 



PACS numbers: 45.70.-n, 61.43.Fs, 65.60.-|-a, 83.80.Fg 



I. INTRODUCTION 



The rapid trend towards the miniaturization of electro- 
mechanical systems has spurred a flurry of theoretical 
activity aimed at investigating quantum mechanical 
effects in the context of more classical subjects like 
heat transfer and mechanical stability (T], 0. Common 
building blocks for these devices are carbon based 
materials, such as nanotubes and graphene, essentially 
two-dimensional elastic sheets which are often bent or 
wrinkled [3, 0] ■ Much effort has been directed towards 
understanding how quantum surface states are affected 
by the underlying curvature of these spaces 0-13 ■ 
However, less attention has been devoted to the converse 
problem: can quantum mechanical effects modify the 
shape and mechanical stability of nano-structures? 

These questions have so far been relegated to the 
fringe of mainstream engineering applications, since 
it is challenging to probe experimentally the regime 
where the characteristic energy of quantum effects is 
comparable to the bending energy. Furthermore, the 
interplay between a quantum order parameter and 
geometry is a more subtle theoretical problem than 
its classical counterpart. Quantum mechanical degrees 
of freedom live in an internal space distinct from the 
local tangent plane of the underlying substrates. As 
a result, coupling mechanisms between the in-plane 
quantum order parameter and curvature are typically 
less intuitive. By contrast, classical buckling is the 
paradigmatic example of how elastic stresses in a crystal 
or liquid-crystal monolayer are screened by curvature 

[IMl- 

In this letter, we demonstrate a quantum analogue 
of buckling and trace its distinctive physical and 
mathematical origins. Specifically, we consider the two 
dimensional (2D) order parameter, tp{p) = |-0|e'^*^''^ 
that describes superfluid or superconducting phases of 
a quantum condensate and show that, the presence of 
isolated point defects (vortices or anti-vortices) where 



the amplitude of the quantum order parameter vanishes, 
causes the substrate to buckle, apparently resembling the 
buckling of flexible membranes by classical crystalline 
defects. However, there is a crucial qualitative difference 
between the two phenomena (summarized pictorially in 
Fig. 1) that stems from the distinct physical mechanisms 
responsible for buckling. Buckling induced by defects 
screening requires that, positive (negative) disclinations 
(eg. 5 or 7-fold coordinated atoms in a crystalline 
monolayer) induce buckling into conical singularities 
of positive (negative) Gaussian curvatures (T9l - [2l| : 
screening requires that the sign of the defect charge is 
matched by the sign of the curvature (see Fig. 1 a-b) . In 
contrast, our stability analysis reveals that despite the 
lack of an explicit defect screening mechanism, quantum 
defects (whose flow lines are shown in Fig. Ic and d 
respectively) will induce buckling into a conical singu- 
larity of negative Gaussian curvature, independently of 
the sign of the defect charge. 

Consider for simplicity, a 2D flexible liquid mem- 
brane of vanishing surface tension supporting a quasi 
two-dimensional quantum condensate. The total energy 
functional H can be written as the sum of contribu- 
tions from the membrane bending energy H,^ and the 
condensate kinetic energy Hy, respectively given by 



(1) 



K 



Here, u = {u^,u^} is a set of two dimensional co- 
ordinates specifying positions R(u) in the plane of 
the membrane, gap = 9qR • and g — det{gai3) 

denote the metric tensor and its determinant, M(u) is 
the extrinsic mean curvature [2^ . k is the membrane 
bending rigidity, K = is the superfluid stiffness 

constant given in terms of the atomic mass m of the 
superfluid and density ps, and da9{u) gives the local 
superfluid velocity. 
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In order to grasp intuitively the origin of the quantum 
buckling instability, consider Eq. © for the case of 
an isolated vortex of topological charge g = ±1 at the 
tip of a conical singularity: an azimuthally symmetric 
surface denoted by a height function (out of plane shift) 
h(p, (j)) — mp, where we have used 2D polar coordinates 
u = {p,(j)}. Note that, the Gaussian curvature for this 
surface is a delta function that vanishes everywhere 
except at the tip of the cone. Due to the azimuthal 
symmetry, we expect the elastic variable = q4> to 
retain its flat space form, so that V0 = ^e^, where is 
the angular unit vector in polar coordinates. With the 
metric expressed in terms of the slant length of the cone 
I — Vl + rri^p, we can evaluate Eq. ([2]) to obtain 



= TrKq^\/l + m? In 



(3) 



where, ao is a microscopic cut-off length (of the order 
of the vortex core radius) and R » ag, is the size 
of the membrane. We see that the energy required 
for the vortex (anti-vortex) to occupy the tip of a 
conical singularity is always greater than its flat-space 



counterpart by a positive definite factor ^/l 



in Eq. 



([3]) 26]. This simple calculation demonstrates that it is 
not energetically favourable for the vortex to buckle the 
membrane into a surface with a positive delta Gaussian 
curvature (the positive definite bending energy in Eq. 
([T]), only adds an extra penalty). 

Contrast the result obtained in Eq. ([3]) with the 
corresponding one for a liquid crystal membrane. In 
the liquid crystal case, the order parameter describes 
the orientation of a vector, not the (scalar) phase of a 
wave-function. This distinction implies that the elastic 
variable must explicitly couple with the underlying 
curvature and thus, each instance of da0 in Eq. ([2]) 
appears in the form da0 — Aa, where the connection Aa 
is a geometric gauge field whose curl equals the Gaussian 
curvature. As a result of this difference, the leading 
order correction in m, appears with a minus sign in 
Eq. and it can be sufficiently large to overcome the 
bending energy cost [l^ H^l . This is the mathematical 
mechanism responsible for the classical buckling of liquid 
crystal as well as crystalline membranes: the geometric 
gauge field couples elastic defects to Gaussian curvature 
via cross terms in {da0 — ^a)^, thereby providing a 
direct mechanism for screening the defect charge. 

Can there be buckling in the absence of a geomet- 
ric gauge field? In order to answer this question, we 
take a more systematic and versatile approach to track 
the curvature correction estimated in Eq. ([3]). For 
general surfaces and defects, curvature corrections arise 
from a conformal anomaly obtained by puncturing the 
domain of integration around defect cores in Eq. ([2]), 
such that the energy remains finite. As a result, the 
vortex free energy becomes the sum of two contributions, 
= Ef + Eg, where Ej is the fiat space energy of the 




FIG. 1: (Color online) Illustrative plots for membrane buck- 
ling corresponding to : (a) positive disclination,(b) negative 
disclination, (c)vortex and (d)anti- vortex. Positive (negative) 
topological defects are shown with red (yellow) dots, where a 
positive (negative) disclination is screened by a surface with 
positive (negative) Gaussian curvature in the top row, while 
a vortex (anti-vortex) buckles the underlying membrane into 
a negative Gaussian curvature surface in the bottom row. 



vortex and. Eg = —Kq'^V{u) is a self-energy contribu- 
tion expressed in terms of the geometric potential V{u.) 
that satisfies a covariant form of Poisson's equation: 



^^^"^(u) = G{u) 



(4) 



Here, the negative of the Gaussian curvature G(u) plays 
a role analogous to the electrostatic charge density @, H ■ 

We now solve Eq. Q for a family of surfaces 
whose height function is described by h{p^(j)) — pf{(j)). 
The geometric potential for a singular distribution of 
Gaussian curvature reads 



V{p) = -27rsF(p + ao,p) 



(5) 



where, the integrated Gaussian curvature s can be ob- 
tained using the Gauss-Bonnet theorem p^ . 
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and the Green's function F evaluated at the core 
of the defect takes the form r(p -|- ao,p) = 
limp^„„ -yl (in ^ - In A) = A In f [la . The coefla- 
cient A is then given by [l9[ 



A = 



2Tr 



Jo [l + /2(0) + /'2(<^)]l/2 



(6) 
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Here, primes denote derivatives with respect to <f). 
Adding the flat space energy Ej — irKq^ In ^ to the 

self energy, Eg ~ —irKq^V = iir'^Kq^Asln'^^, we ob- 
tain the total vortex energy for a singular distribution of 
Gaussian curvature 



E^ ^ TiKq^n + 2ttAs) In ( — 



(7) 



As a check, upon evaluating Eq. ([7]) for a simple 
cone (positive delta Gaussian curvature) described by 

1 



A 



and s — 1 



the result obtained in Eq. ([3|. 



Vl-|-m2 ' 



we recover 



Next, consider a saddle: a surface with negative 
delta Gaussian curvature. We take the simplest surface 
described by a height function, h{p,(j)) — mp cos{2(p). 
The derivation for the total vortex energy proceeds ex- 
actly as outlined above, substituting for f{(j)) = cos(2(/)) 
into the general expressions obtained in Eqs. (HUT]). 
Further, the membrane bending energy can be evaluated 
from Eq. (HI where, for our assumed height function, the 
mean curvature takes the form 12011 . 



M : 



2p[l+/2(0)+/'2(</>)]3/2 



(8) 



Since this surface is not azimuthally symmetric, the resul- 
tant expressions can only be expressed as integrals over 
(j). However, restricting ourselves to small deviation from 
flatness, i.e., m < 1, we can expand 7?^, A and s in a 
perturbation series in m and integrate the resulting ex- 
pressions to obtain the following form, correct to order 
0{m^) for the combined bending and vortex energies 



E = nKq^ 
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(9) 



Here, we have introduced a dimensional parameter, 
r = that serves to quantify the competition between 
condensate kinetic energy and membrane bending energy. 

Inspection of Eq. ([SJ reveals a critical r-c ^ |, be- 
low which the total energy of the buckled membrane 
(m 7^ 0) is less than its flat counterpart (m = 0). 
Further, due to the quadratic dependence on defect 
charge q, this result is independent of the sign of the 
vortex defect. Hence both, vortices and anti-vortices 
will buckle the membrane into a saddle shape. 

In order to determine the shape of the buckled 
membrane, we expand Eq. (|9]) to order 0{m^) and 
evaluate m* corresponding to a minimum in the to- 
tal energy. In Fig. [21 we plot the total energy E 
so obtained, normalized by 'kK\h(^^^, against the 

parameter m for different choices of r. As the ratio is 
decreased below the critical value Tc = 0.17, the energy 
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FIG. 2: (Color online) Normalized energy E{m) — 
E/KttK \n{-^)] as a function of parameter for r = |rc 
(small dots black), r = Tc (solid line black), r = |rc (solid 
red), r — (dash-dot green) and r = (dot blue), where 
Vc ~ 0.17 is the critical ratio of membrane bending rigidity 
K to superfluid stiffness K, below which the membrane can 
buckle. For r < Vc, the normalized energy curves show a min- 
ima below l(top). Bottom plot shows E{m*) as a function of 
r also for the saddle surface, where m* is the value of m for 
which E(m) has a minimum. The plateau for r > rc = 0.17 
shows that the membrane cannot buckle for a ratio higher 
than the critical value. 



(corresponding to red, green and blue curves) has a min- 
imum at TO* indicating that, buckling of the underlying 
membrane into a saddle shape is energetically favourable. 

We next test the validity of our ansatz for the shape of 
the buckled membrane by expressing the height function 
h(p,(j)) — [oi cos(2z0) -|- 5i sin(2i(/))], in the form 

of a truncated Fourier series. Thus, we numerically 
seek the coefficients {ai,bi} that minimizes the total 
Hamiltonians appearing in Eqs. (HJH)), by approximating 
integrals with finite sums. Note, the height function 
used in the perturbation expansion method corresponds 
to retaining just the first term ai in this series. The 
lower panel in Fig. [21 shows the energy so obtained 
as a function of r, where we have used the notation 
E{m*) to denote the energy obtained by solving for 
the Fourier coefficients that minimize the Hamiltonians 
even if we are using more than one term in the Fourier 
series expansion. Above r^, the total energy remains 
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constant, indicating that for a large ratio there is no 
choice of coefficients that minimizes the total energy 
and, therefore, membrane buckling is not favourable. 
However, below Tc, the total energy decreases as r 
is reduced, thus confirming for both a vortex and an 
anti-vortex, the buckling of the underlying membrane 
into a saddle with the dominant contribution to the 
energy coming from only the first term in the Fourier 
series, in agreement with our analytical ansatz. 

Although classical screening of defects by curvature 
of matching sign is absent for the condensate, inspection 
of Eq. (|9]) suggests that the reduced vortex energy, E, 
(dressed by bending contributions) can lower the critical 
temperature, Tc, for Kosterlitz-Thouless transition on a 
flexible substrate. This is apparent from estimatin g T c 
by balancing the vortex energy with its entropy, S, [23|. 
For a two dimensional membrane of size i?, the number 
of possible positions of a vortex is still proportional to 
(^)^ (as in flat space). Any geometric corrections are 
su?)-leading in the limit of a large system size and one 
recovers the familiar result S ~ fc^ In(^)'^, where ks is 
the Boltzmann constant. By contrast, deformation of 
the underlying metric (m ^ 0) changes the pre-factor of 
the logarithmic divergence of the energy E in Eq. 
At the transition, the free energy F = E — TcS = 
vanishes, leading to a reduced critical temperature 

Tc = (10) 

For instance, if the bending rigidity is chosen to be 
below the critical value so that r = ^rc, we can read-off 



from the inset of Fig. ([2]) the corresponding shape of 
the saddle ~ 0.21 yielding, ~ OAlbnK, whereas 
for a flat membrane with m = 0, the corresponding 
temperature will be Tc ^ 0.50ttK, with K of the order 
of 40 Kelvin for a superfluid film of thickness around 
100 Angstroms. 

To conclude, we have demonstrated a striking manifesta- 
tion of the effect of collective excitations (arising from a 
conformal anomaly in the kinetic energy functional of the 
condensate) on the mechanical stability of a membrane 
supporting the condensate. Unlike the prototypical 
buckling instability resulting from the screening of a 
classical defects charge by curvature, quantum buckling 
is a global feature, forced by an energetically more 
favourable target metric. The symmetry of the buckled 
shape with respect to defect charge, corroborates a lack 
of curvature-defect charge screening mechanism and 
is a consequence of the absence of a direct coupling 
between scalar quantum order parameter and geometry. 
In addition, we find that mechanical buckling lowers 
the temperature of the KT transition, providing us 
with an additional experimental observable for quantum 
buckling. Possibilities for future work could include, 
studying the effect of multiple interacting defects and 
the resultant formation of ripples on free standing 
membranes [l^, . 
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